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The present invention relates to an encoding apparatus for two-dimensionaUy 
encoding user data of a user data stteam into channel data stream along a two-dimensional 
chamiel strip of at least two bit rows one dimensionally evolving along a finrt direction and 
being aligned with each other along a second direction, said two directions constituting a 
two-dimensional lattice of bit positions. The present invention relates further to a 
coirespcnding encoding method and a computer program for implementing said method on a 
computer. 

In one-dimensional recording on optical discs the channel bits of the channel 
data stream are recorded along a spiral track, the spiral being one bit wide. For two- 
dimensional recording the channel bits of a channel data stream can also be recorded along a 
spiral, albeit a broad spiral, that consists of a number of bit rows which are aligned with 
respect to each other in the radial direction, that is, in the direction orthogonal to the spiral 
directioa 

European patent application 01203878.2 (ID 607133) discloses a method and 
system for multi-dunensionaUy coding and/or decoding an information to/from a lattice 
structure representing channel bit positions of said coded mformation in at least two 
dimensions. Encoding and/or decoding is perfoimed by usuig a quasi close-packed lattice 
structure. For tiie case of two-dunensional encoding and/or decoding, preferably a quasi- 
hexagonal lattice structure is to be used. 

Further. European patent application 02076665.5 (PHNL 020368) discloses a 
method of multi-dnnensionally encodmg a user data stream of user wonis mto a channel data 
stream of channel words evolving in a one-dimensional direction of mfinite extent Further, a 
corresponding method of decoding is disclosed. The code disclosed in European patent 
application 02076665.5 (PHNL 020368) is based on a finite-state machine. In order to 
implement certain two- or multi-dimensional coding constraints and coding geometries which 
lead to higher storage densities and improve tiie coding efficiency, the metiiod of encoding 
according to tiie fmite-state machine encoder, comprises the following steps: 

a user word is encoded into an NRZ channel word by selecting said NRZ 
channel word from a code table depending on said user word and the cuirent state of tiie 
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underlying finite-state-machine. wherein an NRZ channel word comprises a sequence of 
NRZ channel symbols of NRZ channel bits having a one-dimensional inteipretation along 
said one-dim««ional direction and wherem states of the underlying finite-state-machine 
describing the chamcteristics of the multi-dimensional code are defined by NRZI channel bits 
of the previous chamiel word and by NRZ chamiel symbols of the cuirent chamiel word. 

the NRZ channel symbols are transcoded into NRZI channel symbols by a 
one-dimensional IT-precoding operation including an Integration modulo 2. said IT- 
precoding operation being carried out along said one-dimensional direction of infinite extent 
and 

said finite-state-machme is put into a new state selected fiom said code table 
depending on said user word and the current state of said finite-state-machine together with 
encoding a user word into a channel word. 

In the one-dimensional case a modulation transfer function (MTP) detemmes 
the amount of inter-symbol interference and the achievable Imear density of an encoding 
apparatus along the track direction. The modulation code is designed for data evolving in the 
direction along the track in such a way that the spectral content of the data is transmitted 
properly by the low-pass optical chamiel. Therefore, the MTF is considered as a on^ 
dmiensional function. For the two-dimensional case, however, data evolve in all directions m 
the same way. The shape of MIP is equally important in all directions when a 2D modulation 
code is designed. TTxus. the MTF must be considered as a two-dimensional function, for 
instance as the auto-coirelation functionof a circular top-hat function havmg aparticular 
radius. i.e. where the base of the MTF is a cncle which can be treated as the equivalence of a 
cut-off fiequency beyond which no information is transferred by the chamiel. The 
interpretation of the MTF as the auto-correlation of a circular top-hat that represents the 
detection aperture used m optical recording is described in J. Braat. '«Read-out of Optical 
Discs", in "Principles of Optical Disc Systems", pp. 7-87. Adam Hilger Ltd.. Bristol 1985 • 
This MTF is typical for so-caUed "Central-Aperture" detection, which is preferably used for 
two-dimensional optical storage. 

For two-dimensional high-density storage, it has now been found that the two- 
dimensional frequency space can not be covered fully by the bases of several modulation 
transfer fimctions. i.e. by drawmg a number of touching circles. The two-dimensional 
frequency space consists of a periodic repetition in 2D of the fundamental domain in 
frequency space (which is called the 2D Nyquist cell). For a 2D hexagonal bit-lattice, the 2D 
Nyquist ceU has an hexagonal shape. The circle which is the base of the MTF will at high 
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eaough density not completely cover the full 2D Nyquist cell. As a consequence, there ^ 
fi^uency components in the signal that are not transmitted by the modulation transfer 
function at a combuiation of chosen density and optical pa^meters. Moreover, it has been 
found that th«e are a number of further so-called wonrt case patterns that correspond to these 
frequency components, having a high probability to be detected erroneously when 
transmitted across an optical chamiel for two-dimensional optical storage. 

It is thus an object of the present invention to provide an encoding apparatus 
and a corresponding encoding method which avoid the above described prablem and ensure 
that only data patterns are generated by the encoding apparatus and method which have a 
high probabUity of being detected correcdy by a decoder after being transmitted across an 
optical chamiel. which can. for instance, be a transmission network like the internet or a 
record can-ier like a CD or DVD or any other chamiel with a hard cut-off for the modulation 
transfer characteristic. 

This object is achieved according to the present invention by an encoding 
apparatus as claimed in clarni 1 which comprises a modulation code encoder for modulation 
code encodmg the user data into chamiel data according to a two-dimensional modulation 
code bemg adapted to prevent predetemuned wor^ case patterns of chamiel data in the 
channel data stream. 

The present invention is based on the idea to avoid the generation of such 
worst case patterns as soon and as much as possible during modulation code encoding 
Therefore, these worat case patterns are predetemiined in advance or conditions are defined 
which correspond to the e«stence or generation of such worat case patterns which are known 
to the encoding apparatus so that the encodmg apparatus is equipped with appropriate 
measures so that generation of any wor^ case pattern is prohibited. Preferred embodunents of 
such worst case patterns and appropriate measures for prohibiting the generation hereof are 
defined in the dependent claims. 

According to one aspect of the invention wor^t case pattems shaU be 
prevented that include local periodic patterns havmg fiindamental frequency components 
outside the cirele formed by the base of the two-dimensional optical modulation transfer 
function. Periodic patterns are defined as pattems that have a base on a hexagonal, orthogonal 
or any other lattice, i.e. having base vectors which are linear and integer combinations of the 
base vectors of the hexagonal, orthogonal or any other lattice. These (local) periodic pattems 
are also denoted as super-structures. Upon Fourier transfonnation of these (local) periodic 
pattems a periodic lattice is obtained m Fourier space. Frequency components outside the 
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circle are then understood as points of this Fourier lattice that are outside the circle fonned by 
the cut-off frequency of the modulation transfer function. For obvious reasons of efficientiy 
high data storage, the circle is as small as possible with respect to the extent of the 2D 
fundamental domain . but preferably (m the hexagonal case), has a size such that the circle 
just touches on the fimdamental hexagon (or 2D Nyquist cell) that is the result of the Fourier 
transfonn of the hexagonal lattice itself and is completely contained in this fundamental 
hexagon. Since worst case patterns havmg frequency components outside the circle have a 
high probability of not being detected by a decoder, these shaU be prohibited during 
encoding. 

Among those worst case patterns are patterns which lead to zero HF signal at 
synchronous sampUng moments for a particular target response. Such patterns have circular 
symmetry, for instance a one in the middle and three suirounding one-bits and three 
surroundmg zero-bits. Because the channel is circular symmetric this always leads to the 
same HF signal, and no modulation is present at the output. 

According to another aspect of the invention worst case patterns shall be 
prevented that form a closed ring of bit symbols having alternating bit values or fomi an open 
rmg of bit symbols havmg alternating bit values at a boundary of the chamiel data stream 
where the openmg of the ring is at the side of the boundary of the 2D storage area. These 
patterns are characterized by the feet that the Euclidian weight of the transmitted eiror pattern 
across the optical chamiel is smaU. These kmds of worst case patterns have a high loss with 
respect to the matched filter bound. 

In a preferred embodiment of the invention the modulation code encoder 
comprises a state-transition machine, also called state-transition diagram, describing the flow 
of chamiel words in the direction of evolution of the chamiel code, for prevention of said 
worst case patterns by checkuig the building up of a worst case pattern during modulation 
code encoding and by truncating a sequence of chamiel words building up the start of a worst 
case pattern by entering a state of the encoding apparatus and its state-transition diagram that 
forbids the emission of a NRZ chamiel symbol that leads to a further continuation of the 
worst case pattern. Thus, it is very effectively and veiy soon detected if a worst case pattern 
that shall be suppressed is being built up so that it will not be generated completely, but the 
building up is interrupted. To achieve this, an appropriate state-transition diagram is used. 

A preferred embodiment of such a state-transition diagram is defmed in claim 
5 according to which the state-transition diagram comprises a state conversion unit for 
putting said diagram state-transition diagram into a new state generated depending on said 
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NRZ channel symbol and the cun^t state of said state-transition diagram together with 
encodmg a user word mto a channel word, said generation method being based on an 

of S-m-s^ havmg an identical fen-out except for NRZI channel symbols that lead to a 
5 next NRZ chamiel symbol of a worst case pattern. Additional states are thus added to the 

nonnalSTD-states Of the state transition diagrammonler to generate the extended 
— dxagram on which thecodeconstnu^on will be based, said^^^^^ 
no fen-out for NRZI chamiel symbols that further buUd up worst case patterns 

It is further prefened that the modulation code encoder comprises a chamxel 
10 word conversion unit for transcoding the NRZ chamrel symbols into NRZI chamxel symbols 
by a one-dimensional IT-precoding operation mcluding an integration modulo 2. said IT- 
preceding operation being carried out along said one-dimensional direction of mfinite extent 
Such a channel word conversion unit and a one-dimensional IT-precoding operation are 

15 17^L7^ZZ7 '"'^^^^^^°'°*^^^°---«-'»E-P-Patentapplication 
15 02 076 665.5 (PHNL 020368). 

It can also be advantageous that the modulation code encoder is adapted for 
unplementing a radial and/or tangential k-constraint by adding further STO-states in the state 
transition diagram having no fen-out for predetermined NRZI chamiei symbols. 

A tangential k-constramt is similar to the k-constraint in conventional ID 
10 -^'-^onencoding.Tliek-constraintlunitsthenumberofNRZze^^ 

NRZ ones, which indicate the transitions in the chamiel bitstieam, to a maximmn of k 
T^er^fore. the longest aUowed runlengths ar. Icf 1 bits long. Tto limitation on the occurrence 
of very long r^s is needed since the latter lead to maccuracies in the tuning ..covery. which 
IS dealt with by a device called a phase-locked loop (PLL). The PLL regenerates the internal 
5 clock" that is matched to the length of the bits on the medium: the bit clock is adjusted at 
each occurrence of a transition. Areas on the medium with too few transitions may cause 
"dock-drift". So. m order to retrieve the bit-clock, the receiver needs some changes in the 
bitstream at a large enough pace. For instance, ux the EFMPlus code used in DVD k=10 has 
been used. Shnilar requirements for timing recovery hold in 2D optical storage baled on a 
) bn,ad spiral. In order to retrieve the tangential bit-clock, the number of NRZ 0-symbols 

(..her. one symbol comprises a bit for each bit-row in the strip to be encoded) is limited to a 

^alnumber.For the radial direction, that is orthogonal to the directionofevoluti^^ 
broad spiral, it is advantageous to have a similar k-constramt resulting in changes of the bit- 
patterns in the radial direction. 
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Advantageously the modulation code encoder is adapted for use of 

of «.e codes that .late to prohibiting said wo^t-case patten, is ve^ high. Therefore, the 
code-mappmg from user bits to channel bits typically requi^ the use of ve^ long user 

wordsandve^longchannelwordstoachieveacodewithhighefficienc^^ 

can a^ereforecons^tofalarge number ofNRZ channel sjonbols. Since p„^^^^ 

word xsquxtelargeCfor instance 152 bits forachannel word of 153 channel bits. thusSl 
Channel symbols foras-row^p), it is neariyimpossib^ 

^•'^--'i;-tablein^efonnofamemo,y:thismakestheuseofacon^^^^ 

10 — "-ofcode-tablesdependentontheunderlyingfinite-statemachinecompletely 
impractical^ Therefore, the Channel symbol is generated bas^ 

algonthm that:s called enumerative encoding wWch is la,own per seCseeTJ^.^^^^ 
Enumerative Source Coding". IEEE Transactions on Information Theory, Vol. IT-19 No 1 

15 :™-^^^^^-™------outinre.-t.eandcanbe.^^^^^ 

^'"^^«'"^^'th«^"vention is applicable in a code having a high code rate 
for mstance in a 152 to 153 code for the case If -9 w««, 

annhVH !n *u a However, of course, the invention can be 

applied m other codes as well. 



20 



be explamed in more detail with reference to the 
drawmgs lUustratmg various embodiments of the invention. Therein. 

Fig. 1 shows a block diagram of the general layout of a coding system 
25 scheme ^ ^ '"^^^^^8 ^ ^«p-based two-dimensional 

Fig. 3a shows a schematic diagram of a ID modulation transfer function 
Fig. 3b shows a schematic diagram of a 2D modulation transfer function! 
Fig. 4 shows a linear model for hexagonal distributed bits. 
Fig. 5 shows a first embodiment of a woi^t case pattern, ' 
30 Fig. 6 shows a hexagonal (real-space) bit-lattice with its' base vectors 

^'«-^*°^''*«^iP^°^ lattice with various shapes of possible fundamental 
domams. hexagonal fundamental domain (hfd) and rectangular fimdamental domain (rfd) 

Fig. 8 shows a reciprocal lattice, a fundamental domain and MTF-circles 
representing the 2D cut-off frequencies of the channel. 
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Fig. 9 shows another embodiment of a worst case pattern. 

Fig. 10 shows a table indicating four initial STD-states with NRZI channel 



30 



Fig. 1 1 shows the transition from a current to a next NRZI triplet, 

Fig. 12 shows a table listing some worst case patterns. 

Fig. 13 shows a particular example of a sequence of STD-states. 

Fig. 14a-c shows transitions from a cum«t to a next STD-state in response to 
different NRZ channel symbols, and the elimination of specified symbols in order to truncate 
the generation of a worst-case pattern. 

Fig. 15 shows a first embodiment of a state transition diagram according to the 



invention. 



mvention. 



Fig. 16 shows a connection matrix according to the invention. 

Fig. 17 shows a comiection matrix according to another embodiment of the 



Fig. 18 shows the procedure used during enumerative encoding. 
Fig. 19 shows a flow diagram for enumerative encoding. 

Fig. 20 shows a block diagram of the layout of the coding system in which the 
mvention can be applied. 

Fig. 21 shows a block diagram of an encodmg apparatus according to the 



20 invention. 
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invention, 



and 



Fig. 22 shows a flow diagram followed during enumerative decoding. 
Fig. 23 shows another embodiment of worst case patterns according to the 

Fig. 24 shows a first group of minimum weight patterns, 

Fig. 25 shows a second group of mimmum weight patterns. 

Fig. 26 shows a third group of minimum weight patterns. 

Fig. 27 shows a worst case pattern (a) and its filtered version (b). 

Fig. 28 shows the spectral content of the worst case pattern shovin in Fig. 27a 

Fig. 29 shows an ideal 2d modulation transfer function. 



Fig. 1 shows typical coding and signal processing elements of a data storage 
system. The cycle of user data from input DI to output DO can mclude interleaving 10. error- 
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co^txcn-code (ECC) and modulation encoding 20, 30.signalprep^^^^^^^ 
on fte roooMu^ sO, signal post-processing 60, binary detection 70, and decoding 

80. 90 of the modulation code, and of the interleaved ECC. The ECC encoder 20 adds 
redundancy to the data In order to provide protection against error, from various noise 

. ---^e^C-encodeddataarethenpassedontoamoduladonencoderSOw 

the data to the chamxel. i.e. it manipulates the data mto a form less likely to be corrupted by 
channel errors and more easily detected at the cham^el output. Hre modulated data are then 
input to a recording device, e.g. laser beam recording device or electn,n beam recording 
device, a spatial light modulator or the like, and stort^d in the recordmg medium 50. On the 
retnevmg side, the reading device (e.g. photo-detector device or charge-coupled device 
(CCD)) returns pseudo-analog data values which must be transformed back mto digital data 
(one bit per pixel for binary modulation schemes). The frrst step b this process is a post- 
processing step 60. called equalization, which attempts to undo distortions created in the 
l^corxiing process, still in the pseudo-analog domain. Then the array of pseudo-analog values 
. converted to an arr^y of binary digital data via a bit-detector 70. The array of digital data is 
then passed first to the modulation decoder 80. which performs the inverse operation to 
modulation encoding, and then to an ECC decoder 90. 

In the above mentioned European patent application EP 01203878 2 the 2D 
constrained coding on hexagonal lattices in terms of nearest-neighbour clusters of chamrel . 
bxts rs described. IT^erein. the focus has been mainly on nearest-neighbour constraints with 
therr advantages in terms of mor« robust transmission over the channel, but not on the actual 
construction of such 2D codes. The latter topic is addressed m the European patent 
application 02076665.5 (PHNL 020368). i.e. the unplementation and construction of such a 
2D code is described therein. 

The present invention describes 2D channel coding with different coding 
constraints, the purpose of which is to prohibit certain classes of worst-case patterns. By way 
of example, a certain 2D hexagonal code shall be illustrated in the foUowtag. However it 
should be noted that the general idea of the invention and aU measures can be applied ' 
generally to any 2D code, in particular any 2D hexagonal or square lattice or quasi-closed- 
packed lattice code. Finally, the general idea can also be applied to multi-dimensional codes 
possrbly with isotropic constramts. characterized by a one-dimensional evolution of the code 

Accordmg to the present invention data is recorded on a two-dimensional 
lattice. This lattice is preferably hexagonal, but can also be orthogonal or have any other 
possible periodic repetition of bit positions. The data can be read out by scamring the two- 
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dimensional lattice in a one-dimensional way using an anay of laser spots, preferably with 
the scamiing direction along the direction of the broad spiral. In a preferred embodiment of 
the present invention, the code evolves along the same one-dimensional diiection as the 
scanning direction of the array of laser spots. A broad sphal of the 2D fonnat may consist of 
several 2D strips. Channel code wonls use the full radial extent of a strip. A 2D strip consists 
of a number of ID rows, stacked upon each other in a second direction orthogonal to the first 
du^ction. The principle of strip-based 2D coding is shown in Fig. 2. There can be strips of 
more than one row in the radial direction. Strips within a bn>ad spiral are stacked one upon 
the other, possibly in an alternating fashion when there are two types of strips. Between the 
successive revolutions of the broad spiral, a guard band of. for instance, one row may be 
located. The optical chamiel resulting from the scamiing of the data with each of the laser 
spots has the well-known low-pass behavior with a hard cut-off at f,=2NAA. The 
Modulation Transfer Function (MTF) can be written as: 

m the one-dimensional case this modulation transfer fiinction detennmes the 
amount of inter-symbol mterfeience and the achievable linear density of the system along the 
track direction. In the radial direction the density is detennined by the track-pitch. Tl»e value 
of the track-pitch is determined by optical cross-talk when readmg the data or thennal cross- 
talk (cross-erase) when writing data on a phase-change medium. Hie modulation code is 
designed for data evolving in the direction along the track in such a way that the spectral 
content of the data is transmitted properly by the low-pass optical channel. Therefoi^. only 
the MTF as a one-dimensional function as shown in Fig. 3a is considered. For the two- 
dimensional case however, data evolves in all directions in the same way. The shape of the 
MTP is equaUy important in all dkections when a 2D modulation code is designed. The MTF 
must be considered as a two-dimensional function as is shown in Fig. 3b. The radius of this 
'Chinese haf-function is again ^=2NAA. The two-dimensional MTF can be seen as the auto- 
correlation function of a circular top-hat function with radius NAA. 

In the following the hexagonal lattice will be exemplary considered because it 
IS has a clear advantage in density over a square lattice. The packing fraction of the square 
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lattice is 150/0 smaller than that of the 2D hexagonal lattice. Bits a^ distributed across this 
hexagonal lattice as pits with a diameter smaUer than the hexagonal lattice constant This is 
needed to avoid folding of the signal. The HF signal can be calculated with a scalar 
diffiaction model, but it is also possible to approximate it with reasonable accuracy with a 
hnear model as is shown in Fig. 4 wth a central tap-value equal to co and with c, as tap-value 
for all 6 nearest-neighbour bits around the central bit 

At a first glance the most critical (highest) frequency on this lattice seems to 
be the periodic repetition of zeroes (-) and ones (+) as is mdicated in Fig. 5. From the Imear 
model a HF signal is obtained equal to: 

HF = ±(Co~2c,) 

where the sign depends on the polarity of the cluster. 

The HF signal is equal to zero if co=2c, for all bit positions of the critical 
pattern. This mdicates that at the cut-off of the MTF the linear model fulfills this condition of 
co=2c.. The hexagonal lattice parameter aH determines the distance between the planes of 
zeroes and ones for the critical pattern. It is equal to d=V3 aH. The frequency resulting from 
tiiis period d has a cut-off at 2NAA: 

20 1 = — L_=2J5t4^ __1 x_ 

Inserting X =405 nm and NA=0.85 leads to a lattice panmieter of 137.5 mn It 
can thus be expected that a scalar diffi^ction sunulation with this lattice parameter leads to 
HF levels that reasonably fit the linear model for which co=2c,. This appears to be true. Hie 
25 equivalent areal density is then 2.2 times the areal densitiy achieved in BD (and 2.0 times the 
density of BD after inserting one so called guard-band for every II tracks in the broad spiral 
each guard band has the width of one bit row). 

A second check for the above calculation is derived from the calculation of 
reciprocal lattices. To do this the hexagonal lattice is described by two base vector: 



30 
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This is shoMm schematically in Fig. 6. The equivalent of such a vector in the 
ID case is a scalar that denotes the lattice constant (spatially) or the sampling period T (for 
time). 

Now. a so-called generating matrix can be defined (for a giveb 2D periodic 
5 lattice, L, tiiis generating matrix is not unique): 



10 



15 



20 



25 



[ei, 62,} 2 [ 1 _ij 



The generating matrix of tiie reciprocal lattice is the inveise transposed of the 
generating matrix of the origmal lattice (the lattice upon which the bits a« o«iered; also 
known as real-space or du-ect-space lattice) : 



Oh iS 



The reciprocal lattice is also a hexagonal lattice. It is rotated by nl6 witii 
respect to tiie original lattice and it is scaled by a factor of I/2V3 as is schematically shown in 
Fig. 7 (It IS to be noted tiiat only the hexagonal shape of the lattice is rotated by nl6 The 
individual base vectors are rotated by nl6 and -nie i^spectively). TT,e figure also shows some 
possible fimdamental domains of the reciprocal lattice. A fundamental domain is an area ti.at 
has the property that adjacent fimdamental domains do not overlap and tiaat ti.e concatenation 
of fundamental domains of all fl»e lattice points must fill tiie complete 2D space. The 2D 
Shamion samplmg theorem now states tiiat a function in tiie spatial domain can be retrieved 
from its sampled vei^ion if tiie Fourier Transform is contained witiiin a fundamental domain 
of the r^iprocal lattice. There is no unique fundamental domam as is indicated m Fig 7 but 
preferably tiie hexagonal fundamental domain is chosen (it is no coincidence tiiat such a 
hexagonal domain is tiie best fit to tiie circular shape of tiie 2D MTF. which supports tiie 
above argumentation m terms of packmg fractions). It shaU be noted further, tiiat tiie 
rectangular fundamental domains in Fig. 7 show tiiat for a hexagonal lattice a somewhat 
higher frequency can be allowed m one direction if tiie frequency in tiie orthogonal direction 
30 IS known to be limited to a somewhat lower value. 

The distance between points on tiie reciprocal lattice is 2/V3aH. This means 
tiiat ,f touchmg circles around each lattice point can be drawn tiiey have a radius of l/VSaH. 
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The ID equivalent of this radius can be treated as a sort of cut-off frequency of the channel 
The circles can be seen as the base-plane of the MTF that is shown in Fig. 3b. The radius of 
such a circle must be equal to 2NA/X. This means that it can be written: 

The fact that this value is equal to the value that was aheady derived above 
verifies the reasoning. Besides the lattice itself it is also possible to define generation 
matrices of periodic patterns on this lattice (these are know as super-structures m 
crystallography). The critical pattern discussed above for example can be written as: 



LeiU 



^£»^^/3 Vs] 
2[l -3j 

The corresponding reciprocal lattice is: 
2 1 p 1 1 



The vectors of this lattice are drawn in the reciprocal lattice. This is shown m 
Fig. 8. Furthermore, the circles as discussed above are drawn. It now appears that indeed one 
of the reciprocal vectors of the critical pattern Ues exactly on the circle (i.e. exactly on the 
cut-off of the MTF). This is indicated with a small cross (in the duection of fa). 

In the above touching cucles in the frequency domain have been drawn that 
can be seen as the base of the MTF. These circles were treated as the equivalence of the cut- 
off frequency. It appears that the frequency space cannot be fiiUy covered by drawing 
touchmg circles. This means that there might be frequency components m the signal that are 
not transmitted by the modulation transfer function at the combmation of chosen density and 
optical parameters. When investigating this problem it indeed appears that there are bit 
patterns that have most or all of their energy in the area between the circles. The most 
promment example is the one mdicated in Fig. 9. The base vectors of the lattice are also 
shown in Fig. 9. It has reciprocal vectors on the comers of the hexagonal fundamental 
domain of the reciprocal lattice mdicated with a circled cross in Fig. 9. 
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The HF-samples from this pattern are equal to co-6c, or -co. This means that 
the HF-samples do not show any modulation when gomg from one bit of the critical pattern 
to the other when: 

co-6ci = -Co „ Co = 3ci 
This shows that the chamiel is already 'dead' for lower densities than the one mitially 
expected. This also means that if the worst-case pattern occure for a longer time the PLL of 
the encoder might fell out of lock or shows severe phase deviation. 

One solution to the problem is to lower the density m such a way that the 
circles overlap and cover the complete fi^quency space. However, this will cost quite some 
density. The increase in radius of the circle must be a fector l/cos(30°)=2/V3 as is clear from 
Fig. 8. The density ratio is simply calculated as {2l<3f =4/3. TOs means that the achievable 
density is only 75% of the initially assumed density i.e. 0.75*2.2=1.65 times the densitiy of 
BD. This is, however, not yet a sufficient density. 

A second, more desirable possibility is to exclude patterns from the data that 
have frequency components in the non-covered region. This can be done by applying a 
suitable modulation code. A known modulation code is described in European patent 
appHcation 02076665.5 (PHNL 020368). Therein a 2D modulation code is proposed that 
evolves in one dunension along a strip consisting of rows (a so-called fish bone code). 
Each of the channel words consist of a sequence of M-aiy (M=2»'««') NRZ channel symbols 
that indicate the type of transition between the successive NRZI bits that are finaUy recorded 
on the medium. This modulation code was described by considering the practical case of 3 
rows in a 2D strip. For simplicity it will also be adhered hei« to this assumption although the 
present mvention is generic for any number of rows larger than 3. Tlie concept of isolated 
bits that was introduced in the above mentioned European patent application is not needed, 
and therefore only 4 states are required in the initial state transition diagram (STD). These 
states are listed in Fig. 10. 

Since it is not necessary to distmguish between the different polarities of a 
(worst-case) pattern the states are denoted with x and y indicating opposite binary symbols 
(i.e. when x is defined as 0 or 1, y is defmed as 1 or 0 respectively). The reason that the exact 
polarity must not necessarily be considered is that the output data, which are preferably 
recorded on a disc, is of the bipolar NRZI type which is generated from the NRZ data by 
usmg a 1 -T precoder (integrator modulo 2). In the present case the 1-T precoder operates in a 
1-dimensional way along the direction of the spiral. As already mentioned, a state transition 
can be described by a chamiel symbol, in the present case consistmg of a bit-triplet where 
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each of the NRZ bits mdicate a transition (1) or the absence of a transition (0) in the bipolar 
NRZI channel bit-stream. This is shown in Fig. II. 

Now it is possible to write down the sequence of NRZ channel symbols that 
wiU generate the worst-case-pattem that was shown in Fig. 9. It must be noted that there are 
other wor^-case-pattems. Some of them are listed (not exhaustively) in the table shown in 
Fig. 12. 

Long patterns can be formed by repeating these basic buildmg blocks. TTiey 
will not show modulation for the central row for a particular choice for the ratio co/c,. When 
started in state 03 the sequence for the pattern of Fig. 9 is as shown in Fig. 13. 

It should be noted that the last state is equal to the first state and the sequence 
can be repeated to form a longer worst-case pattern. It is desm«i to eUminate the sequence of 
NRZ symbols that lead to the worst-case pattern by truncating this sequence as soon as a 
possible build up of the pattern is detected. In the followmg example the sequence is aheady 
truncated after two consecutive NRZI triplets that are part of the worst-case pattern (of course 
also longer sequences might be allowed depending on the robustness of the bit-detection for 
this pattern and the allowed burst-like behaviour of random errors in the system). The 
tnmcation can be done by defining 'colored' states that have the same NRZI triplet but do not 
have a fan-out of a symbol that will lead to the next triplet of the worst-case pattern. This is 
shown schematically in Figs. 14a to 14c for difiTerent initial states. 

In the above cases 05 is a version of as without fen-out for the NRZ symbol 
(5). m the same way is a version of a, without symbol (2) and a, is a version of 03 without 
symbol (7). This means lhat the state transition diagram (STD) now looks as shown m Fig. 
15. The corresponding connection matruc is shown in Fig. 16. 

The largest eigenvalue. of this connection matrix is equal to 7.91697. This 
25 means that the capacity of this code is equal to: 

C=ilog2Un«)= 0.994983 
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The fector 1/3 is needed because each state corresponds to a channel symbol that consist of 3 
bits. It can be seen that the code rate is nearly equal to 1 . This is not smprising because the 
occurrence of the worst-case pattern was not expected very frequently. As already mentioned 
above an additional problem of the worst-case pattern is that the modulation is zero and the 
PLL or delay synchronization block (for the aligmnent of the signal waveforms obtamed 
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from the array of laser-spots that is inclined with respect to the direction of the broad spiral) 
might experience problems. Other causes for problems in the PLL or delay synchn,nization 
block are a long number of consecutive zeroes or ones. In view of robust timing recovery it 
IS beneficial to include a k-constramt A k-constmint can be applied in both tangential K, Ind 
5 radial direction K,. To introduce a tangential k-constraint the maximum number of 
consecutive symbols (0) that occur between non-zero symbols can be limited to K.. For 
example when K,=l the encoder is allowed to emit a single symbol (0). Then it will be 
arrived in an extra state a, From this state it is not aUowed to emit a symbol (0). Such a 
construction means that the number of states in the state transition diagram doubles This 
10 results in the comiection matrix shown in Fig. 17. The same procedure can be appUed for K. 
constraints larger than 1 . 

It can be seen that on the main diagonal of the states o, to 04 the value is 
reduced by 1 in comparison to above comiection matrix shown in Fig. 16. TTiis is caused by 
the feet that emitting a symbol (0) in these states will not lead to the same state, but will lead 
15 to an extra state a,' to o,'. This causes I4 (the 4x4 identity matrix) in this part of the 

comiection matrix. When a symbol (0) is emitted in the states a, to a, this will automatically 
stop the generation of a worst-case pattern. As a result, the next state will be one out of o. ' to 
C4'. As an example, when symbol (0) is emitted in os the next state will not be os' but a,' 
When being in the states o,' to a,' any symbol can be emitted again except the symbol (0) 
20 Therefore, this part of the matrix is a copy of the upper left comer of the same matrix and the 
rest is fiUed up with ze«,es. For K.=2 another 7 states can be added and the same procedure 
can be foUowed. These are the code constraints for which the practical code construction is 
described in the following deswiption. 

The radial constramt K. might be violated in state a, because here all bits of 
25 the triplet are equal. A transition to a subsequent state a, can be via a symbol (0) or a 
symbol(7). The symbol (0) was already excluded by the introduction of a tangential 
constraint K.. Therefore implementing the constraint K. can simply be done by mtroducing an 
additional state that is equal to a, but has no fan-out for the symbol (7). TTie radial k- 
constraint might become of more importance when also re-sampling of the data in the radial 
30 direction (in case of sloppy radial trackmg) is considered. 

Because the capacity of the chamiel constraints with K,=2 is nearly equal to 1 
in a preferred embodunent a 152 153 code mapping is chosen. This results in 51 fishbones 
of3 bits. 1 fishbone is added forDC-control on one ofthe bit-rows in the strip of3 bit-rows 
(comparable to BD where 1 DC-control bit is added every 45 source bits). The 152 source 
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bits lead to 152/8=19 bytes. When the code is implemented that eliminates the woi^t-case 
patterns it appears that indeed the worst-case patterns do not occur within the chosen strip of 
3 rows. However, concatenation of the strip might again lead to the worst-case pattern (e g 
when strip A and B are concatenated a woist-case pattern is still possible in a 3-row strip that 
IS buUd up from the bottom row of strip A and the two top rows of strip B). Therefore, a 
separation track (which can also be considered as a strip of one bit row wide) is prefer^ly 
inserted in between the strips that are 3 bit-row wide. The code for this separation track has a 
k-constraint and eliminates both the 'lOOlOOlOOl...' and '01 101 101 10...' sequences (NRZI) 
that might lead to the worst case pattern. This can be done by applying a 12 -* 13 code (the 
code construction is not discussed here, but the well-known ACH procedure with state- 
splitting and state-mergmg can be followed). 4*13 perfectly fits the number of fishbones 
(51+1) and perfectly fits a byte oriented approach: Each user word of 12 bits is 1.5 bytes and 
4*1.5=6 bytes. So this example has atotal of 69 (3*19+2*6) bytes encoded m one block with 
1 1 rows. Another number of rows (equal to 4n+3 with n an integer number) is also possible 
with the same fonnat, but it has to be taken into account that in the present proposal the 
central row has to be DC-free for servo purposes (and the 12 - 13 code is not DC-free): this 
means the current 12-> 13 code is applicable for the cases where n is even (3, 1 1, 19,... bit- 
rows). 

The implementation of the 152 153 code normally leads to unpractically 
la^ge look-up tables (or code-books). A possible way to cu-cumvent this problem is the 
application of enumerative coding as described in T.M.Cover "Enumerative Somce Coding- 
IEEE Trans. Inf Theory. Vol IT-19. No. 1. Jan. 1973. Each possible codeword is given an ' 
index according to its position in the lexicographic ordering (with lexicographic ordering it is 
simply meant that for example in the decimal system a 4 has higher order than 0. 1. 2 and 3 

andmthebinary systema 1 has higher order than 0 or in other notation: 0<1). For a binary 
sequence the index can be calculated as: 



30 



. . n 

H-^J= X xjnJx},X2, xj-),0] 

where ns(xi, x^, . . Xk) is the number of sequences for which the first k coordinates are given 

"^^^ ^"^Ple shown in Fig. 1 8 illustrates the procedure. In the code 
implementation for the current invention, binary sequences are not obtamed, but sequences of 
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M-aiy symbols (in the present embodiment M=2'*-»=8). The symbol is an element of 
2, 3, .... M-1}. The lexicographic mdex can now be written as: 
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y = 1 m = 0 

The sequence of M-axy symbols can be considered as the chamiel data and the 
lexicographic index of each of the symbols in the sequence is the equivalent modulation 
decoded sequence. Code tables are not required. 

For the encoding the inverse procedure can be used. It is started with data that 
indicate the lexicographic index. It shall be assumed that the lexicographic index is(x) = 1 is 
given. THe following steps must be followed for the binary case (i.e. the inverse of the above 
equation for is(x)): 

Stepl: Set X to zero; 

Step2: If i>ns(0) set xi=l and set i=i-ns(0); 

Step3: For k=2. .... n. if i>ns(x,. X2, .... x^.,. 0) set Xk=l and set i=i-ns(x,. X2. 

....Xfc.i, 0). 

Following this procedure will result in the encoding (in case of binary data) A 
similar encoding procedure can be followed for M-aiy data in the 2D case. Also here no large 
code tables are needed. According to the above last equation for is(x). only the set of 
numbers denoted by ns(x,, x^, . . xj.,. m) is needed. In the d-constrained ID RLL binaiy 
coding, this set of numbers reduces to xj Nd(n-j) (with xj the binaiy bit value, with n the 
channel-word length, and with Na(p) the number of d-constrained sequences of length p). A 
similar reduction cannot be made for the current case of the 2D coding in a strip of 3 bit-rows 
aligned on an hexagonal lattice of bits. 

Supposing that encoding is started from state So always being equal to a, of 
the state-transition diagram (the state at moment j is indicated with 5^; this state is equal to 
one of the states in the state transition diagram. i.e. S,- e {o„ c„ .... a„}). What is actually 
needed is to take mto account at symbol position j for emission of symbol m is the resulting 
state, hereby considering the histoiy of the previously emitted chamiel symbols x,. x^ x , 
which represent a path through the state-transition diagram. This is shown schematic^ly ' 
below, with the resulting state denoted by Sj": 
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X.^E.-^22-JeU£3....JizlUx>.^2:;. 

It is thus obtained 

n^xi,xi xy- !,«)=«, 

5 

vAere n's^ indicates the total fen-out ftom state Sj" up to any possible state q after 
completion of the channel word (with emission of n-j extra symbolsuntil the end). Therefore 
it is obtained that 

where D is the comiection matrix of the state-transition diagram, the required power is n-j 
(the remainmg number of symbols to be emitted), the required matrix element is (Sj™. q) and 
the sum is over all possible states q. 
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The algorithm for enumerative encoding is shown in more detail in Fig 19 
With the user-wotd is associated a unique index. TTie user-wonl is mapped onto a chamiel 
word consisting of n fish-bones. With each fish-bone coiwsponds an M-aiy chamiel symbol 
(with M =2**-). For the practical case that is considered here, it is obtained that N„v^3, and 
M=8. The M-aty chamiel symbols a=l. 2. .... n) for the successive fish-bones in the chamiel- 
word are computed one by one^ In the ve,y fi„rt step (SI), the mitial state So (prior to the first 
symbol to be generated) is set to the first state al of the state-transition-diagram. 

Then, a loop is started for all possible symbols m, startmg fit)m 0 (S2). With 
each symbol xD] is also associated a flag denoted by flag[j]. which is initially set to zero at 
the start of the m-loop. Next (S3), it is checked whether candidate symbol m for symbol 
position j can be emitted from state 2,.. and whether the flag of the current symbol position 
IS not yet raised to "1 If this is true ("Y"). then a temporary parameter "itmp" is computed 
(S4) which is the number of possible chamiel sequences that start with symbol m at position 
J. For that puipose, first the next-state is computed which results from emission of candidate 
symbol m at position j: this state is referred to as E,-. The number of remaining symbols to be 
emitted in order to reach the end of the chamiel word equals n-j. and the total fan-out from 
state Lj» to any possible state q at the end of the chamiel word is computed as a sum over q of 
matrix-element (Sj". q) of the (n^).th power of tiie comiection matrix D (S5). This fan-out 
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depends on three parameter: (i) the state to start from, thus taking into account the history of 
the previous symbols; (ii) the position] of the channel symbol under investigation; and (iii) 
the candidate symbol m that is considered as channel symbol for position j. If the index (or 
what is left of it after previous iterations, denoted as "jtmp") is larger than the parameter 
"itmp" (S6). then "jtmp" is decreased by "itmp" (S7), and the symbol at position j is (at least) 
equal to m+1; in the other case ("jtmp" being smaller than "itmp") the flag for the cunent 
position j is raised to "1" (S8). which impUes that no further action is done on this position 
for larger values of m. 

After completion of the m-loop (S9), it is checked whether the resulting 
symbol xD] equals "0" (SIO), and m such case, whether this is allowed &om state Sj.i in view 
of the Kt constraint (SI 1); if not, then the symbol is raised by one, thus xD]=l (S12). This 
yields the final decision for xQ]. Then the finally resulting next-state 2j is computed (S13). 
Next the position "j" is increased by one (S14). and if not yet the end of the channel word has 
been reached ("Y", thus j < n) (S15), then it is going back to start with the m-loop for the 
new current position j. In case j=n+I ("N"). the end of the channel word has been reached, 
and all channel symbols have been determined (denoted by (x[l], x[2], x[n])). This 
completes the encoding procedure. 

Another possible implementation for the elimination of worst-case patterns is 
the use of -forbidden- sites. With forbidden sites those bit positions are meant that will lead in 
combination with previous data bits to a worst case pattern when they are occupied with a 
particular bit-value, independent of the actual bit-value that was planned to be written at that 
location. These forbidden bit-positions are just filled with a bit-value that prohibits formation 
of the worst-case pattern, and the next bit-position in the 2D channel stream is considered for 
wilting the next bit. Due to the need to skip some forbidden bit-positions, there is a certain 
expansion area needed to accommodate storage space for the skipped positions. Similar to 
forbidding the fen-out of a particular symbol a bit-position is marked as -forbidden'. Because 
it is not possible to store data on these sites, the capacity of the channel would become of a 
stochastic nature, depending on the number of worst-case patterns that occur. This can be 
solved by adding a small part for reserved data, which is the expansion area. The procedure is 
now as follows: NRZI bits are written in a sequential way through the (for instance 3-row) 
strip; when a certain bit-position would lead to a worst-case pattern for 0 or 1, it is written 
with the bit-value that does not generate the worst case pattern. The data that can not be 
stored because of these forbidden sites is added in the reserved data block. When no 
forbidden sites are used at aU. the reserved expansion area is filled with zeroes. The size of 
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the reserved data block equals the maximum number of forbidden bits. The data (including 
bits in the reserved part) is used as input to the error correction block. Parities are calculated 
that also contain information about the bits in the reserved block. The data (without the 
reserved block) and parities are now transmitted across the channel. At the receiver side the 
data are not avaUable that was in the reserved part These data can be restored by setting it 
initially to zero and using the ECC at cost of efficiency of the ECC (2 parity symbols for each 
•missing' byte of reserved data). In the (veiy) rare case that the ECC is not able to conect a 
block and to retrieve the reserved data because it needs too many parities for the reserved 
block and to correct random errors, the ECC processing can be repeated by settmg erasures 
(which requires 1 parity for each 'missing' byte) in the reserved block instead of zeroes 
(which will cost 2 parities for each 'missing' byte). The complete procedure is schematicaUy 
shown in Fig. 20. 

A more detailed block diagram of an encoding apparatus is shown in Fig. 21. 
The whole system runs on a clock mclk. First a user word (user index) Uk is loaded and the 
counters j and m are reset to one and zero respectively. The initial state So is reset to 1. Based 
on this initial state and a first value of m a next state is calculated S j". Based on this state a 
value itmp is calculated (see the flow diagram in Fig. 1 9) in the icalc block. This value itmp 
is subtracted from the value jtmp (which is initiaUy equal to U^). The result Otmp-itmp) is 
checked if it is negative or not. If not a next value of m is evaluated (triggered by mclk). Tliis 
again gives a new value of itmp. This is agam subtracted from jtmp Otmp is now the updated 
value equal to jtmp-itmp of the previous iteration). If now the sign of Gtmp-itmp) is smaller 
than zero the value of jtmp is taken via mux 2 as the input to muxl to prevent that itmp is 
subtracted from jtmp. In feet, jtmp-itmp was only used as an intennediate value to check the 
sign. Now the output of the Sign determination block is 1 and the following wiU happen: 
^- is set to m. This is one symbol of the channel word. 

2. Mcounter is reset to zero for the determination of the next channel symbol 

xD+1]. 

3. Jcounter is incremented for determination of the next channel symbol xD+1]. 
As soon as the jcounter is larger than n a new user word is loaded and the jcounter is reset 
The final result (a series of xQ] values for j=l to n) is the output in the form of the channel 
word. It should be noted that this block diagram does not include the k-constraints for a 
limited number of zero symbols. 

A flowchart of a decodmg method is shown in Fig. 22. The userword given by 
the user mdex is first set to zero (S20). Then the channel word (code word) is taken as an 
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input. The first symbol of the codeword is taken (j=l) and an index m is iterated from 0 to 
xQ] (S21). Each time for a new value of m the next state is calculated (S22) and from this 
state the total number of possible sequences that start with symbol m at position j are 
calculated (S23). The user index is mcremented with this number (S24) and then the index m 
5 is incremented by 1 (S25) until it reaches a value larger than x\jl At that moment (S26) the 
next symbol from the channel word is taken G=j+1). If it was the last channel word (S27) 
then the calculated index is the final user word. 

Above different worst case patterns for the 2D optical channel have been 
described. A worst case pattern is a pattern that has a high probability to be detected 
1 0 erroneously when transmitted across the channel. For example when a=l is transmitted and 
a'=-l is received the error is -2. In the following the error sequence is divided by 2 to obtain a 
more compact notation: e = (a - a')/2. Further, a more detailed overview and a catalog of 
these patterns shall be given. For the time being two different types of worst case patterns 
shall be distinguished: 

Type-I: Error patterns that have frequency components beyond the cut-off of 
the 2D-MTF. These patterns have been described above in detail. 

Type-II: Minimum weight error patterns. 
It shall be noted that the definitions are not orthogonal (e.g. the minimum weight error 
patterns have then" main spectral content at high frequencies and even part of it outside the 
2D-MTF). The classification is based on the way the patterns were found. 

Type-II error patterns are characterized by the fact that the Euclidian weight of 
the transmitted error pattern across the chaxmel is small: 
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where g is the channel response. These kind of patterns have a high loss with respect to the 
matched filter bound. When a brute force search is done it has been found that all worst case 
patterns have the form of a closed rmg of alternated +1 and -1 error symbols. A graphical 
representation is given in Fig. 23. The stars indicate a don*t care symbol. At the boundaries 
30 the ring might be open at the side of the boundary. A catalog of these type of patterns is 
shown in Figs. 24 to 26. 

To eliminate the possible occurrence of these kinds of error patterns in the 
channel bitstream the states that are the result of these error patterns and/or the states that 
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allow these error patterns to occur must be eliininated.Code constraints are needed in such a 
way tiiat an occurrence of these types of errors in an admitted data sequence (i.e. one that 
does not violate the coding constraints) leads to a violation in the coding constraints. When it 
is looked carefully at the pattern it is found that the only patterns that allovwj the error 
5 patterns of Fig. 23 are the inverted form of the patterns itself. Therefore, it is sufficient to 
eliminate all sequences that form a closed ring of alternated +1 and -1 symbols fiom the data 
pattern. Fortunately, for the worst case one indicated in the upper left comer of Fig. 23 this is 
automatically done vnth the code that was described above. For larger rings (that have a 
smaller loss with respect to die matched filter bound) this is not the case. The fact that these 
10 patterns have such a high loss with respect to the matched filter bound also becomes clear 
from the hexagonal Fourier transfonn of these patterns. The example of Fig. 27 results m the 
spectrum of Fig. 28. 

Looking at this spectrum and comparing it to the transfer function of the 
channel that is given in Fig. 29 it can be immediately seen that most of the spectral content of 
15 the pattern is in a range where the transfer of the channel is only small. A considerable part of 
the spectral content is even beyond the cut-off of the optical channel. 

In conclusion 2 ^es of patterns can be distinguished that must be eliminated 
fi-om the data sequence (based on the way the patterns were found). These can be defined as 
follows: 

local periodic patterns (superstructures) that upon 2D Fourier transformation 
of the lattice formed by tiie pattern result in a lattice in the Fourier space with points beyond 
tile cut-off of the modulation transfer fimction of tiie optical channel. 

patterns tiiat form a closed rmg of alternated +1 and -1 symbols or part of a 
ring at tiie boundary of a stiipe, where tiie opening of tiie ring is at tiie side of tiie boundary of 
25 the broad spiral. 
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